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ABSTRACT 

This  research  has  two  main  components:  (1)  a  new, 
wavelet-based  cascade  model  for  intermittent  (uneven)  dis¬ 
tributions  of  objects  in  terrestrial  environments  and  (2)  a 
complementary  analysis  of  intermittent  wave  scattering. 
Taken  together,  these  developments  facilitate  realistic  sta¬ 
tistical  descriptions  of  the  performance  range  (uncertainty) 
associated  with  sensing  and  communication  systems.  The 
wavelet  cascade  model  uses  a  simple  set  of  decay  processes 
to  generate  complex  random  media  with  realistic  organiza¬ 
tion  of  objects  in  space  and  time.  It  is  used  to  recreate  the 
complex  geology  of  the  lava  flow  field  around  the  Amboy 
Crater,  CA,  and  incorporated  into  simulations  of  seismic 
wave  propagation  for  this  site.  The  simulations  reproduce 
well  the  signal  travel  time  distributions  recorded  along  a 
sampling  of  864  separate  propagation  paths. 

1.  INTRODUCTION 

Terrestrial  environments,  both  natural  and  man-made, 
often  possess  complex  spatial  distributions  of  objects  that 
are  intermittent  {uneven,  irregular ,  bunched ,  or  sporadic). 
The  distribution  of  volcanic  rock  in  Fig.  1,  from  the  Am¬ 
boy  Crater,  CA  vicinity,  is  an  example.  Besides  surface 
geology,  intermittency  is  exhibited  by  atmospheric  turbu¬ 
lence,  weather  events,  vegetation,  animal  populations,  ur¬ 
ban  construction,  and  many  other  natural  and  man-made 
environments.  This  paper  endeavors  to  provide  a  consis¬ 
tent,  general  approach  to  describing  intermittency  and  its 
effects  on  an  important  Army  problem:  signal  transmission 
through  terrestrial  environments.  The  performance  of  sens¬ 
ing  and  communication  systems  is  heavily  impacted  by  the 
environments  in  which  they  operate.  Reliance  upon  such 
systems  will  continue  to  grow  with  advancement  in  infor¬ 
mation  collection,  distribution,  and  analysis  technologies. 

Conventional  treatments  of  acoustic,  seismic,  and  elec¬ 
tromagnetic  (EM)  signal  transmission  have  assumed  that 
heterogeneities  (terrestrial  objects)  responsible  for  random 
scattering  of  the  wavefield  are  evenly  distributed  in  space. 


This  is  a  significant  shortcoming  in  existing  capabilities 
for  characterizing  and  predicting  environmental  effects  on 
system  performance.  Figure  2  depicts  how  an  intermittent 
medium  affects  transmission  and  sensing  of  waves.  Along 
some  signal  propagation  paths,  the  waves  encounter  many 
scattering  objects  and  signal  coherence  (a  measure  of  sig¬ 
nal  randomness,  which  usually  determines  the  sensing  sys¬ 
tem  performance)  is  strongly  degraded.  Other  paths  en¬ 
counter  few  scattering  objects  and  thus  provide  high  co¬ 
herence.  Conventional  wave  scattering  theory  (Tatarskii, 
1971)  assumes,  in  effect,  that  all  paths  experience  similar 
scattering.  It  thus  neglects  the  range  of  possible  outcomes 
and  related  uncertainty  in  system  performance,  which,  de¬ 
pending  on  the  environment,  can  be  very  significant. 

The  research  described  here  involves  two  main  com¬ 
ponents:  first,  a  new,  wavelet-based  cascade  process  model 
for  describing  intermittent  distributions  of  objects  in  ran¬ 
dom  media,  and  second,  a  complementary  description  of 
the  effect  of  the  intermittency  on  signal  transmission.  Al¬ 
though  the  discussion  emphasizes  surface  geology  and  seis¬ 
mic  wave  propagation,  the  goal  is  to  formulate  a  concep¬ 
tual  and  statistical  model  for  intermittency  and  its  effects 
on  wave  propagation  with  applicability  extending  well  be¬ 
yond  the  geologic  domain.  In  particular,  the  model  is  also 
intended  for  application  to  atmospheric  turbulence,  which, 
by  scattering  sound,  radio,  and  optical  waves,  often  plays 
a  key  role  in  the  performance  of  sensing  and  communi¬ 
cation  systems.  Intermittent  signal  propagation  in  urban 
areas  is  also  a  consideration:  acoustic  or  optical  sensors 
may  be  rendered  useless  by  obstructions  in  built-up  sec¬ 
tions,  whereas  they  may  produce  very  clean  signals  around 
open  areas  such  as  large  lawns  and  parking  lots. 

This  paper  is  structured  as  follows.  In  Sec.  2,  we  de¬ 
scribe  the  underlying  conceptual  model  for  intermittency 
and  its  appearance  through  cascade  processes ,  which  in¬ 
volve  a  sequence  of  self-similar  reactions  by  which  one 
generation  of  objects  is  created  from  the  next.  Section  3 
describes  a  basis  for  mathematical  construction  for  random 
media  models,  called  quasi-wavelets  (QWs),  that  is  partic- 
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Fig.  1:  Intermittent  distribution  of  volcanic  (basaltic) 
rock  in  the  Amboy  Crater  vicinity,  Mojave  Desert,  CA.  The 
crater  itself  is  visible  on  the  horizon,  in  the  middle  of  the 
picture. 

ularly  well  suited  to  describing  intermittency  and  cascade 
processes.  Section  4  brings  the  paper  back  to  the  inter¬ 
mittent  surface  geology  shown  in  Fig.  1:  a  QW  model  for 
the  complex,  random  geology  at  this  site  is  described,  and 
simulations  of  seismic  wave  propagation  through  this  sta¬ 
tistical  model  are  compared  to  seismic  observations  from 
a  unique  experiment  that  sampled  wave  propagation  along 
864  paths. 

2.  INTERMITTENT  CASCADE  PROCESSES 

A  cascade  process  begins  with  an  initial  (or  root )  gen¬ 
eration  of  objects  triggering  a  reaction  by  which  the  next 
generation  is  formed.  This  generation  then  triggers  forma¬ 
tion  of  the  next,  and  so  on.  The  cascade  may  involve  ei¬ 
ther  relatively  large  objects  breaking  down  (decaying)  into 
smaller  ones,  or  relatively  small  objects  consolidating  into 
larger  ones.  If  the  cascade  is  self  similar,  the  transformation 
from  one  generation  to  the  next  proceeds  by  application  of 
a  scale-independent  set  of  rules.  To  formulate  such  scale- 
independent  rules,  dimensional  quantities  (such  as  sizes  of 
the  objects,  the  duration  of  the  generation,  and  energy)  are 
normalized  in  some  manner. 

Figure  3  illustrates,  in  two  dimensions,  a  self-similar 
cascade  mechanism  used  by  Frisch  et  al.  (1978)  to  describe 
turbulence.  The  largest  eddies  (those  produced  directly  by 
the  instability  responsible  for  the  turbulence)  are  space¬ 
filling  as  represented  by  the  initial,  solid  square.  Each  time 
the  generational  process  of  decay  into  smaller  eddies  is  re¬ 
peated,  some  portion  of  space  randomly  becomes  inactive 
(unoccupied  by  the  smaller  eddies).  In  this  particular  illus¬ 
tration,  1/4  of  the  space  becomes  inactive  upon  each  iter¬ 
ation.  This  procedure  leads  to  a  fractal  set  (Mandelbrot, 
1977;  Schroeder,  1991)  upon  which  the  dissipation  of  tur¬ 
bulent  kinetic  energy  into  heat  takes  place  (Frisch  et  al., 
1978;  Meneveau  and  Sreenivasan,  1987). 

Although  the  volcanic  processes  originally  responsible 
for  the  formation  of  the  Amboy  Crater  landscape,  and  sub¬ 
sequent  wind  and  water  erosion,  are  quite  different  from 
turbulence,  the  distribution  of  basaltic  material  in  Fig.  1 
bears  a  striking  resemblance  to  the  later  stages  of  the  ran- 


FlG.  2:  Interference  with  propagating  waves  and  sensor 
performance  resulting  from  an  intermittent  material  distri¬ 
bution. 
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Fig.  3:  Notional  cascade  mechanism  used  by  Frisch  et 
al.  to  describe  intermittent  dissipation  in  turbulence.  Each 
panel  indicates  a  generation  of  the  turbulence.  The  cascade 
proceeds  counter-clockwise  from  the  upper  left.  Light  blue 
shading  indicates  active  regions  of  the  turbulence.  Repiti- 
tion  of  the  cascade  leads  to  a  fractal  set. 

dom  cascade  process  depicted  in  Fig.  3.  One  might  hy¬ 
pothesize  that  self-similar  processes  have  been  or  are  at 
work  at  Amboy;  for  example,  lava  ducts  may  tend  to  or¬ 
ganize  in  a  self-similar  fashion  (with  larger  ducts  feeding 
into  smaller  “fingers”),  and  erosion  breaks  down  basaltic 
rocks  in  a  cascade-like  fashion.  Still,  we  do  not  assert  that 
a  turbulence  cascade  model  like  in  Fig.  3  applies  literally  to 
past  or  on-going  processes  at  the  Amboy  site.  Rather,  we 
take  the  similarity  in  appearance  as  evidence  that  models 
for  intermittent  cascade  processes  might  be  suitable  for  de¬ 
scribing  the  geometry  of  other  random,  intermittent  distri¬ 
butions  of  objects  (such  as  urban  construction,  vegetation, 
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clouds,  and  the  many  other  terrestrial  and  galactic  phenom¬ 
ena  described  in  Mandelbrot  (1977)),  even  when  the  link 
between  the  cascade  process  and  intermittent  geometry  is 
not  readily  apparent. 

Let  us  thus  make  several  observations  about  the  cas¬ 
cade  process  depicted  in  Fig.  3.  First,  there  are  certain 
obvious,  unrealistic  aspects  of  Fig.  3  in  comparison  to  the 
Amboy  photograph,  Fig.  1.  (This  is  unsurprising,  since  the 
original  purpose  of  Frisch  et  al.  (1978)  was  to  explain 
the  formation  of  a  fractal  dissipation  surface  in  turbulence, 
rather  than  to  create  visually  realistic  scenes.)  One  obvious 
idealization  is  the  square  shape  of  the  objects  in  the  cas¬ 
cade.  This  may  be  a  problem  in  wave  propagation  model¬ 
ing,  for  example,  when  parallel  and  perpendicular  surfaces 
combine  to  form  artificial  waveguides  and  resonance  cham¬ 
bers.  Note  also  in  Fig.  3  that  the  children  objects  fit  together 
snugly  without  any  overlapping,  into  the  space  occupied  by 
the  parent  object.  Such  a  fit  is  most  easily  accomplished 
with  squares  and  cubes,  and  requires  that  the  number  of 
children  is  an  integer  multiple  of  the  number  of  parents.  In 
many  cascade  processes,  including  turbulence,  the  children 
do  not  possess  sharp  boundaries  and  thus  must  overlap  in 
some  sense;  they  also  have  an  apparently  random  spatial 
distribution.  The  objects  usually  possess  a  continuum  of 
sizes,  rather  than  an  integer  progression.  Figure  3  also  im¬ 
plies  perfect  alignment  in  a  temporal  sense:  all  reactions 
creating  a  particular  object  size  are  synchronous. 

Another  observation  concerns  the  dynamic  nature  of 
the  cascades,  and  how  they  overlap  in  space  and  time.  Over 
time,  erosion  will  break  down  the  rock  at  the  Amboy  site, 
creating  smaller  objects.  Similarly,  turbulent  eddies  decay 
and  are  dissipated  unless  renewed  by  an  instability  mecha¬ 
nism  such  as  shear.  Thus,  the  age  of  the  cascade  is  an  im¬ 
portant  consideration.  The  cascade  may  be  younger  (with 
relatively  larger  objects)  in  some  regions,  and  older  (with 
relatively  small  objects)  in  other  regions.  Different  parts 
of  the  cascade  may  also  proceed  at  different  rates:  smaller 
objects  usually  break  down  more  quickly  than  larger  ones. 
This  implies  the  existence  of  two  mechanisms  by  which  the 
cascade  objects  can  become  progressively  less  space-filling 
as  their  size  is  decreased.  The  first,  illustrated  directly  in 
Fig.  3,  is  that  parents  do  not  have  enough  children  to  oc¬ 
cupy  their  original  space.  The  second  is  that  subsequent 
generations  are  represented  by  fewer  objects,  by  volume, 
because  they  do  not  survive  as  long. 

With  these  considerations  in  mind,  we  now  de¬ 
scribe  a  more  general  framework  for  intermittent  cascade 
processes.  To  start,  we  define  a  number  of  discrete  size 
classes  for  the  objects.  The  characteristic  size  (such  as  the 
radius  or  side  length  of  a  square)  for  each  class  n  will  be 
designated  as  an .  A  smaller  index  represents  larger  objects. 
The  ratio  a  =  an+\  /an  of  one  size  class  to  the  preceding  is 
allowed  to  take  on  any  value,  although  it  must  be  indepen¬ 
dent  of  n  for  a  self- similar  construction. 

Next,  we  formulate  a  set  of  conservative  reactions  by 
which  a  single  parent  object,  size  an ,  breaks  down  into  mul¬ 
tiple,  smaller  objects  with  sizes  an+\,  an+2, _ Each  reac¬ 


tion  is  defined  by  the  number  and  energies  of  the  children 
objects  it  produces.  One  illustrative  reaction  is 

an  ->  (N29pi)an+2.  (1) 

The  notation  here  means  that  an  decays  into  N2  children 
objects  of  size  an+ 2,  each  with  energy  p2  times  the  energy 
of  the  parent  object.  Conservation  of  energy  requires  that 
N2P2  =  1.  The  particular  case  N2  =  3  and  a  =  \/\fl 
corresponds  to  the  reaction  in  Fig.  3.  When  the  objects 
consist  of  a  different  material  than  the  background  medium, 
mass-conservation  constraints  should  also  be  implemented. 
Other  reactions  might  be  defined  that  produce  multiple 
sizes  of  children  objects,  such  as 

an  ->  (N\,pi)an+\  +  (Ni,p3)  an+s.  (2) 

This  reaction  can  produce  a  scale  densification;  that  is,  the 
distribution  of  object  sizes  involved  in  the  cascade  becomes 
finer. 

How  many  reactions  do  we  need,  and  how  realistic 
must  they  be,  to  describe  a  particular,  intermittent  geom¬ 
etry?  Some  perspective  can  be  gained  from  recent  develop¬ 
ments  in  agent-based  modeling  (Bonabeau,  2002),  which 
show  that  complex  group  behavior  can  emerge  from  a  num¬ 
ber  of  autonomous  entities  obeying  simple  rules.  Thus,  it 
might  be  suggested  that  just  a  few,  carefully  chosen  reac¬ 
tions,  for  example  reactions  that  lead  to  cascades  with  ap¬ 
propriate  fractal  dimensions  and  scale  densification,  will 
suffice.  This  idea  is  demonstrated  with  an  example  in 
Sec.  3. 

For  self-similar  dynamics,  the  ratio  of  the  lifetime  zn+\ 
of  one  size  class  to  the  preceding  must  be  constant;  that  is, 
zn+\ /xn  =  r,  where  r  is  independent  of  n.  The  lifetime  can 
be  interpreted  as  the  characteristic  decay  time  for  comple¬ 
tion  of  a  reaction  involving  a  parent  with  size  an .  Various 
probabilistic  models  for  the  decay  may  be  considered.  For 
present  purposes,  we  adopt  the  following  simple  descrip¬ 
tion.  After  creation,  each  object  of  size  an  persists  for  a 
time  zn .  One  of  the  R  defined  reactions  then  takes  place 
with  assigned  probability  pr,  r  =  1, . . . ,  R.  The  pr  must 
sum  to  1 . 

Our  description  of  the  reaction  is  incomplete  with¬ 
out  a  procedure  for  distributing  the  objects  in  space.  This 
procedure  actually  generates  the  intermittency,  if  it  exists. 
Mahrt  (1989)  has  proposed  classifying  intermittency  as  in¬ 
trinsic  or  global.  Global  intermittency  is  associated  with 
the  production  mechanism,  which  in  our  cascade  would  be 
the  distribution  of  the  largest  (root)  generation  objects.  In¬ 
trinsic  intermittency  is  associated  with  the  number  density 
of  the  children  objects  relative  to  the  parents.  The  process 
illustrated  in  Fig.  3  is  an  example  of  intrinsic  intermittency. 
For  the  volcanic  flow  shown  in  Fig.  1(a),  the  large  patches 
of  rock  are  associated  with  global  intermittency  of  the  lava 
flow,  whereas  cooling,  erosion,  and  other  processes  subse¬ 
quent  to  the  volcanic  activity  tend  to  produce  intrinsic  in¬ 
termittency.  In  the  next  section,  we  will  discuss  a  statistical 
approach  for  placing  the  objects  that  includes  both  intrinsic 
and  global  intermittency. 
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3.  QUASI- WAVELET  MODELING  OF 
INTERMITTENCY 

3.1  What  are  Quasi- Wavelets? 

Until  this  point,  the  objects  in  the  cascades  have  not 
been  explicitly  described.  We  now  associate  them  with 
wavelet,  or  wavelet-like,  functions.  Wavelets  make  conve¬ 
nient  analogs  for  physical  objects  because  they  are  spatially 
localized  and  their  mathematical  properties  are  well  under¬ 
stood  (e.g.,  Torrence  and  Compo,  1998).  We  use  the  term 
quasi-wavelets  (QWs)  because,  strictly  speaking,  our  ob¬ 
jects  are  not  always  true  wavelets.  In  particular,  restrictions 
such  as  a  zero  mean  may  be  relaxed,  and  the  QWs  may 
be  spherically  symmetric  in  multiple  dimensions,  which  is 
not  normal  practice  in  wavelet  analysis.  Like  customary 
wavelets,  QWs  are  always  based  on  translations  and  dila¬ 
tions  of  a  parent  function;  however,  their  positions  and  ori¬ 
entations  are  randomly  generated. 

Mathematically,  the  QW  is  defined  by  its  parent 
function  /  (<f),  which  we  assume  here  to  be  spatially 
symmetric.  Here  f  is  the  magnitude  of  the  vector 
£  =  (r  —  bnm )  / an ,  r  the  spatial  coordinate,  and  bnm  the 
center  position  of  the  QW.  The  index  n  indicates  the  size 
class  of  the  QW  and  m  a  particular  QW  within  that  size 
class.  An  individual  QW  comprising  a  scalar  field  (i.e., 
a  single  field  such  as  temperature)  is  given  by  Tnm  (r)  = 
hnmATnf  (|r  _  bnm  |  /an)  , where  hnm  is  a  random  sign 
factor  and  A  Tn  is  an  amplitude  factor.  The  total  field  is 
found  by  summing  the  contributions  of  all  individual  QWs 
comprising  a  size  class  and  then  summing  over  the  size 
classes. 

We  have  considered  a  number  of  QW  parent  functions 
for  various  applications.  Examples  are  shown  in  Fig.  4. 
Three  of  the  illustrated  QWs  are  used  to  represent  scalar 
fields.  To  represent  a  turbulent  velocity  field,  the  parent 
function  can  be  used  as  a  vector  potential ;  the  velocity  field 
associated  with  the  QW  is  determined  by  taking  the  curl  of 
the  potential.  Each  QW  is  thus  analogous  to  a  turbulent 
eddy.  When  the  size  distribution  and  rotation  rates  of  the 
QWs  are  chosen  in  a  manner  consistent  with  Kolmogorov’s 
hypotheses,  turbulence  fields  with  realistic  spectral  proper¬ 
ties  result  (Goedecke  et  al.,  2004).  The  velocity  QWs  can 
be  paired  off  with  scalar,  dipole  QWs  to  represent  corre¬ 
lated,  random  fields.  In  Wilson  et  al.  (2008a),  this  concept 
was  used  to  represent  the  turbulent  heat  transport  above  a 
heated  boundary. 

3.2  Construction  of  Intermittent  Fields  from  QWs 

Previously,  random  fields  have  often  been  constructed 
from  harmonic  functions.  The  familiar  technique  involves 
scaling  the  amplitudes  of  the  harmonic  functions  by  the 
square  root  of  the  power  spectrum,  randomizing  the  phases, 
and  then  taking  the  inverse  Fourier  transform  (Shinozuka 
and  Deodatis,  1996).  However,  such  an  approach  leads  to 
random  fields  with  smoothly  distributed  statistical  proper¬ 
ties:  it  is  not  suitable  for  describing  intermittency.  This  is 


Fig.  4:  Images  of  various  3D  quasi-wavelets  (QWs).  A 
section  of  each  QW  has  been  cut  away  to  show  the  inte¬ 
rior  structure.  Top  left:  Monopole  (spherically  symmetric) 
“top-hat”  QW.  Top  right:  Monopole  Gaussian  QW.  Bottom 
left:  Dipole  derivative-of-Gaussian  QW.  Bottom  right:  Ro¬ 
tating  (velocity)  Gaussian  QW. 


our  main  motivation  for  constructing  random  fields  from 
QWs  instead. 

The  basic  approach  was  described  in  Sec.  2.  We  pro¬ 
pose  here  a  simple,  2D  implementation,  appropriate  to  the 
Amboy  site  shown  in  Fig.  1,  based  on  the  following  three 
reactions:  an  (4,  l/4)aw+2 ,  an  ->  (2,2/5)an+2  + 
(1,  1/5)  an+2,  and  an  ->  (1,  \/2)an+\  +  (4,  1/8 )aw+3. 
Here,  the  ratio  a  =  an+\/an  is  set  to  1/a/2.  The  first  of 
these  reactions  might  be  considered  the  “baseline”  reaction, 
which  involves  decay  of  an  object  into  four,  equal-sized, 
equal-energy  objects.  The  second  produces  a  spatial  fractal 
distribution  and  diversification  of  energy.  The  third  pro¬ 
duces  diversification  of  spatial  scale.  Furthermore,  we  set 
the  lifetimes  to  rn  =  (4/5)n~l  z\ ,  so  that  smaller  objects 
decay  more  rapidly.  (Note  that  the  actual  value  of  z\  is  im¬ 
material.)  At  time  zn,  one  of  the  reactions  r  will  take  place 
with  equal  probability  ( pr  =  1/3).  The  objects  are  repre¬ 
sented  as  solid  circles  (the  2D  version  of  the  QW  in  Fig.  4, 
upper  left).  A  recursive  Matlab  program  was  written  to 
simulate  the  cascade  development. 

Figure  5  shows  a  random  realization  of  this  cascade 
process  based  on  a  uniformly  random  spatial  placement 
rule  for  the  children  objects;  that  is,  the  children  are  ran¬ 
domly  placed  anywhere  in  the  domain,  thus  eliminating  the 
spatial  organization.  A  total  of  16  root  objects,  each  with 
radius  30  m,  were  placed  with  uniform,  random  distribu¬ 
tion  about  the  250  m2  simulation  domain.  (An  edge  region, 
with  width  25  m,  is  not  shown  in  the  figures.)  The  initia¬ 
tion  times  (start  of  the  root  generation)  for  the  16  cascades 
varied  randomly  between  —z\  and  +z\.  The  random  field 
is  shown  at  time  z\  =  2,  i.e.,  when  most  of  the  cascades  are 
mature. 
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Fig.  5:  A  random  medium  constructed  from  objects 
(quasi-wavelets)  evenly  distributed  in  space.  In  analogy  to 
the  Amboy  Crater  site,  black  represents  the  basaltic  rocks 
(cascade  objects);  gold  is  the  sand  (background  medium). 

Figure  6  was  created  in  the  same  manner,  except  that 
the  children  objects  are  placed  near  to  their  parents.  The 
center  coordinates  of  the  children  were  selected  according 
to  a  2D  Gaussian  distribution,  with  standard  deviation  2 an . 
Hence,  the  children  fall  within  a  radius  2 an  of  the  parent 
about  68%  of  the  time.  The  cascade  now  has  pronounced 
intermittency,  in  both  intrinsic  and  global  senses.  There  is 
also  clustering  of  objects  based  on  their  size.  This  occurs 
because  of  the  different  ages  of  the  cascades  dominating  in 
different  regions. 

3.3  Wave  Scattering  by  QWs 

We  have  previously  described  a  theory  for  wave  scat¬ 
tering  by  QWs  (Wilson  et  al.,  2004).  The  basic  idea  is 
to  calculate  scattering  from  a  single  QW,  and  then,  assum¬ 
ing  a  single-scattering  approximation  is  valid,  the  scattered 
fields  from  all  QWs  in  all  size  classes  can  be  summed  to 
obtain  the  total  scattering  cross  section.  Recently  Wilson 
et  al.  (2008b)  used  QW  models  as  the  basis  for  numeri¬ 
cal  simulations  and  analytical  formulas  of  wave  propaga¬ 
tion  through  turbulence  with  intrinsic  and  global  intermit¬ 
tency.  The  results  demonstrate  that  global  intermittency 
leads  to  highly  non-Gaussian  signal  statistics  for  phase 
fluctuations  of  acoustic  signals.  The  non-Gaussian  statis¬ 
tics  have  important  implications  regarding  applicability  of 
the  conventional  theory  of  wave  scattering  (  e.g.,  Tatarskii 
(1971),  which  is  based  on  Gaussian  signal  statistics)  as  well 
as  practical  strategies  for  employing  sensors.  Practically 
speaking,  signal  processing  strategies  should  be  based  on 
the  recognition  that  performance  of  sensing  systems  con¬ 
sists  of  periods  of  good  performance  disrupted  by  episodes 
of  dramatically  degraded  performance.  This  suggests  that 
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Fig.  6:  Same  as  Fig.  5,  except  that  the  objects  are  dis¬ 
tributed  in  space  according  to  a  cascade  process  including 
intrinsic  intermittency. 

“consensus-based”  signal  processing,  as  opposed  to  con¬ 
ventional  cross-correlation  processing,  is  most  appropriate. 

4.  GEOLOGIC  INTERMITTENCY  AND  SEISMIC 
WAVE  PROPAGATION 

4.1  Amboy  Crater  Site 

As  mentioned  in  the  Introduction,  we  conducted  a  seis¬ 
mic  propagation  experiment  in  the  vicinity  of  the  Amboy 
Crater,  which  is  located  in  the  Mojave  Desert  near  Amboy, 
CA,  about  75  miles  east  of  Barstow.  This  site  was  cho¬ 
sen  for  its  highly  inhomogeneous  geology.  The  volcanic 
flow  field  covers  an  area  of  about  70  km2  and  consists  pri¬ 
marily  of  alkali  basalt  that  erupted  onto  the  alluvial  plain 
about  74  to  85  ka  (Phillips,  2003).  The  flow  field  is  com¬ 
posed  primarily  of  vesicular  pahoehoe,  with  a  hummocky 
surface  and  relief  ranging  from  2  to  5  m  (Parker,  1963). 
Aeolian  (windblown)  sand  fills  the  depressions  and  forms  a 
discontinuous  cover  over  the  flow  field.  Most  areas  of  sand 
lack  vegetation  and  commonly  have  irregular,  pebble-  to 
cobble-size  gravel  scattered  across  their  surface.  Sand  de¬ 
posits  generally  range  in  thickness  from  a  few  cm  to  more 
than  2  m.  The  sand  at  depth  greater  than  ~30  cm  is  often 
denser  and  wet. 

The  experimental  seismic  grid  (to  be  described  in  more 
detail  shortly)  was  established  in  a  200  m2  region  of  the 
northwest  section  of  the  lava  flow  field  shown  in  the  Fig.  7 
satellite  image.  A  high-resolution  digital  elevation  map 
(DEM)  and  material  classification  were  developed  for  the 
area  using  a  tripod-mounted,  scanning  lidar,  Riegl  Model 
LMS-Z420i.  The  final  DEM  (see  Fig.  8)  was  generated 
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Fig.  7:  High-resolution  satellite  (QuickBird)  image  of 
the  area  containing  the  seismic  test  grid.  Seismic  sensors 
and  shot  points  are  shown  in  yellow. 


Time  (s) 


Fig.  9:  Signals  recorded  at  the  geophones  approximately 
60-m  distant  from  each  of  the  source  locations. 


Fig.  8:  Digital  elevation  map  for  the  area  containing  the 
seismic  grid.  Open  circles  are  the  geophone  locations;  tri¬ 
angles  are  the  weight-drop  locations. 

with  1-m  grid  spacing.  An  ISODATA  unsupervised  clas¬ 
sification  scheme  used  the  intensity  of  the  lidar  returns  to 
map  the  material  types. 

4.2  Seismic  Measurement  Procedure 

Thirty-six  vertical  component  geophones  (Mark  Prod¬ 
ucts  L-28  with  4. 5 -Hz  resonance)  were  installed  in  a  6-by- 
6  square  grid  with  30-m  spacing.  The  geophone  locations 
were  adjusted  as  necessary  so  they  were  always  inserted 
into  sand.  The  seismic-wave  source  was  a  16-kg  steel 
sphere  dropped  from  a  height  of  1.5  m.  Twenty-four  dif¬ 
ferent  source  locations  around  the  edge  of  the  sensor  grid 


were  used  as  shown  in  Fig.  8.  A  geophone  located  0.5  m 
from  the  impact  location  triggered  the  recording  system,  a 
Geometries  NZ  seismograph.  The  digital  recordings  had  a 
24-bit  dynamic  range  and  a  sampling  rate  of  1  kHz. 

The  36  geophones  and  24  source  locations  yield  864 
propagation  paths.  Figure  9  shows  2-s  time  series  from  a 
single  weight  drop  at  each  source  location,  as  recorded  by 
the  geophone  located  (approximately)  60  meters  from  the 
source.  These  measurements  show  truly  remarkable  dif¬ 
ferences  for  such  a  short  propagation  distance  and  slight 
location  changes.  There  are  large  variations  in  the  shape 
and  arrival  time  of  the  initial  pulse,  as  well  as  the  envelope 
of  the  entire  time  series  and  the  character  of  the  later  ar¬ 
rivals.  In  some  cases,  the  maximum  amplitude  is  early  in 
the  time  series  (0.1  s)  while  in  others  it  is  delayed  up  to 
0.3  s.  The  high-frequency  arrival  at  around  0.17  s  in  some 
of  the  traces  is  likely  the  acoustic  “thump”  of  the  weight 
impacting  the  ground. 

4.3  Travel-Time  Analysis 

We  now  analyze  the  arrival  time  distribution  of  the  sig¬ 
nal  energy.  Figure  9  shows  that  the  signal  energy  arrives 
at  different  times  along  various  propagation  paths,  and  is 
spread  over  a  long  time  interval.  To  quantify  these  effects, 
we  process  the  signals  by  first  taking  the  Hilbert  transform. 
The  magnitude  of  the  Hilbert  transform  describes  the  enve¬ 
lope  of  the  signal,  the  square  of  which  may  be  interpreted 
as  the  signal  energy.  The  signal  energy  calculated  in  this 
manner  is  then  integrated  from  the  start  to  end  of  each  2-s 
record.  The  linear  trend  resulting  from  background  noise 
was  subtracted  from  the  integrated  energy. 

The  results  of  this  analysis  are  shown  in  Fig.  10.  The 
figure  is  a  scatter  plot  of  the  10%,  50%,  and  90%  energy  ar¬ 
rival  times  as  a  function  of  the  distance  between  the  source 
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Fig.  10:  Scatter  plot  of  signal  travel  times  vs.  distance. 
The  travel  times  are  normalized  for  a  propagation  speed  of 
600  m/s  as  explained  in  the  text.  Blue  dots  are  the  times 
at  which  10%  of  the  signal  energy  has  arrived,  green  dots 
50%,  and  red  dots  90%. 

point  and  geophone.  The  times  are  normalized  by  dividing 
by  the  distance,  and  multiplying  by  600  m/s.  This  normal¬ 
ization  “flattens”  the  dimensionless  time  for  the  10%  ar¬ 
rival  time  at  a  value  of  approximately  1,  thus  indicating  that 
the  early-arriving  energy  in  the  recordings  travels  at  about 
600  m/s.  The  normalized  plot  also  makes  clear  that  the 
spread  (duration)  of  the  signal  energy  is  relatively  larger 
for  shorter  paths. 

Due  to  the  complexity  of  the  site  geology  and  seis¬ 
mic  wave  propagation  in  general,  interpretation  of  the  sig¬ 
nals  is  not  immediately  clear.  Since  P-waves  (body  com- 
pressional  waves)  travel  through  basalt  at  5400-6400  m/s 
(Press,  1966),  it  is  unlikely  that  most  of  the  energy  reaching 
the  geophones  propagates  in  this  form.  The  P-wave  speed 
in  unsaturated  sand  varies  widely  over  a  range  of  about 
200-1000  m/s  (Press,  1966).  S-wave  (body  shear  wave) 
speeds  are  typically  60%  of  P-wave  speeds;  Rayleigh  (sur¬ 
face)  waves  in  a  half  space  travel  somewhat  more  slowly 
than  the  shear  body  waves.  Given  the  observed  propagation 
speed  of  600  m/s  for  the  early  arriving  energy,  it  is  reason¬ 
able  to  associate  this  energy  with  waves  traveling  primarily 
through  sand.  Since  the  sand  consists  of  a  thin  layer  with 
variable  thickness,  the  waves  likely  interact  with  both  the 
upper  and  lower  surfaces  of  the  sand  and  thus  are  neither 
purely  body  nor  Rayleigh  waves.  The  typically  long  delays 
between  the  10%  and  90%  energy  arrival  times  suggest  that 
the  waves  are  propagating  along  paths  through  the  sand  and 
reflecting  from  the  basaltic  features. 

4.4  Comparison  to  QW  Scattering  Model 

To  simulate  wave  propagation  at  the  Amboy  site, 
we  use  an  acoustic  finite-difference,  time-domain  (FDTD) 
code.  For  simplicity,  the  FDTD  code  is  acoustic  in 
the  sense  of  propagating  a  single,  compressional  wave. 


x  (m) 

Fig.  11:  Finite-difference,  time-domain  calculation  of 
waves  propagating  through  a  random  model  for  the  Am¬ 
boy  site.  Red  is  positive  pressure;  green  is  negative.  (The 
amplitude  is  arbitrary.) 

Second-order,  central  finite  differencing  is  used  in  both 
space  and  time.  The  source  is  represented  by  a  Ricker 
wavelet  with  center  frequency  50  Hz.  The  sandy  areas 
are  modeled  with  a  phase  speed  of  600  m/s  and  density 
1600  kg/m3.  The  basalt  is  modeled  with  a  phase  speed  of 
6000  m/s  and  density  1900  kg/m3.  Seismic  attenuation  is 
included  in  the  calculations  by  a  recursive  filter  method. 

Figure  1 1  shows  the  FDTD  calculation  for  the  partic¬ 
ular  random  realization  of  the  Amboy  site  shown  in  Fig.  6. 
Shown  is  the  pressure  at  0. 1717  s  after  the  onset  of  the  ini¬ 
tial  source  signal.  Consistent  with  the  seismic  traces  (Fig. 
9),  there  is  very  strong,  random  scattering  of  the  signal.  The 
initial  signal,  which  at  600  m/s  travels  about  100  m  dur¬ 
ing  this  time  interval,  is  strongly  dispersed  and  scattered. 
The  longer  wavelength  within  the  basalt  is  evident  by  the 
smoother  field  within  the  large  pieces  of  basalt. 

Normalized  travel  times  (10%,  50%,  and  90%  energy 
arrival  times,  normalized  for  a  propagation  speed  of  600 
m/s  as  described  earlier)  were  calculated  from  the  simu¬ 
lations  by  placing  receivers  at  5-m  intervals  along  radial 
spokes  at  16  different  azimuths  form  the  source.  Proba¬ 
bility  density  functions  (pdfs)  for  the  travel  times  from  the 
simulation  are  compared  to  the  experimental  data  Fig.  12. 
The  pdfs  are  based  on  combining  all  distances  between  50 
m  and  70  m.  The  mode  (pdf  maximum)  and  spread  of  the 
arrival  times  are  reproduced  quite  well  by  the  simulation. 
Interestingly,  the  mode  for  the  50%  arrival  time  is  some¬ 
what  too  low,  and  that  for  the  90%  arrival  time  somewhat 
too  high.  This  may  indicate  that  the  wave  speed  and  atten¬ 
uation  need  some  adjustment.  We  are  presently  analyzing 
in  more  detail  ground-truth  data  from  the  site,  such  as  the 
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Fig.  12:  Probability  density  functions  for  the  normalized 
travel  times  from  the  Amboy  data  (solid  lines)  and  FDTD 
simulation  (dashed  lines).  Distances  between  50  m  and  70 
m  are  included.  Colors  indicate  the  percentage  of  the  total 
signal  energy. 

lidar  material  classifications  and  sand  depth  measurements, 
so  as  to  refine  the  details  of  the  QW  model. 

5.  CONCLUSION 

Random  scattering  objects  in  terrestrial  environments, 
such  as  turbulence,  geologic  heterogeneities,  vegetation, 
and  man-made  structures,  often  impede  the  performance 
of  Army  sensing  and  communication  systems.  These  ob¬ 
jects  typically  are  distributed  intermittently  (unevenly)  in 
space,  thus  causing  system  performance  to  vary  widely 
from  one  location  to  another.  By  combining  fractal  cas¬ 
cade  processes  and  wavelets,  we  have  devised  a  method 
for  constructing  random  fields  with  spatial  intermittency. 
The  model  distributes  the  originating  (root)  objects  of  the 
cascades  through  space  and  time.  This  is  termed  global  in¬ 
termittency.  The  cascade  process  concentrates  smaller  ob¬ 
jects  into  progressively  smaller  regions  of  space.  This  is 
intrinsic  intermittency.  Using  a  cascade  model  based  on  a 
few,  simple  reactions,  we  were  able  to  create  random  me¬ 
dia  with  a  qualitative  appearance  very  similar  to  a  site  with 
highly  complex  geology,  namely  the  Amboy  Crater  flow 
field.  We  then  simulated  the  statistical  properties  of  travel 
times  of  waves  propagating  through  this  complex  environ¬ 
ment.  Good  agreement  was  found  between  the  observed 
and  simulated  travel-time  distributions. 
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